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Abstract 

We study a Markov process on a system of interlacing particles. At large times the particles fill a 
domain that depends on a parameter e > 0. The domain has two cusps, one pointing up and one pointing 
down. In the limit e J. the cusps touch, thus forming a tacnode. The main result of the paper is a 
derivation of the local correlation kernel around the tacnode in the transition regime e J, 0. We also prove 
that the local process interpolates between the Pearcey process and the GUE minor process. 

Keywords: Determinantal point processes. Random growth. GUE minor process. Pearcey process. 

1 Introduction 

In a recent paper [5] the authors introduced a Markov process on a system of interlacing particles. This 
model contains many parameters, creating a rich pool of interesting particular examples, and at the same 
time it has an integrable structure that allows for explicit computations. It is the purpose of this paper, to 
study a special case of this model. The interest of this model lies in the fact that for large time the particles 
will fill a domain that has a tacnode on the boundary. A somewhat similar situation also occurs in the case 
of non-intersecting Brownian paths with multiple sources and sinks [1]. The local process at the tacnode 
in this model is not understood, although a conjecture is given in pQ. The integrability of the model we 
consider allows us to compute the local process around the tacnode. This is the main result of this paper. 
We consider an evolution on particles that are placed on the grid 

f to + 1 1 

G = \ (x,m) | m = 1,2,... xeZ+— — \. (1.1) 



Hence, if (x, to) <E G, then x takes integer values for odd values of to and half-integer values for even values 
of m. At each horizontal m-section we put to particles and denote their horizontal coordinates by x™ for 
k = 1, . . . , m. The evolution is such that at each time the system of particles satisfies the interlacing condition 

Z/P-i < ^k-x < x ™, k = 2,...,m, to = 2,... (1.2) 

At time t — we put the particles at positions x™ = — (to + l)/2 + k as shown in Figure [TJ 

The evolution of the particles is as follows: each particle has two independent exponential clocks, a left 
and a right clock respectively. If the right (left) clock rings, the particle attempts to jump to the right (left) 
by one. But in doing so it is forced to respect the interlacing condition according to the following two rules: 
if the right (respectively left) clock of the particle at x™ rings, then 

1. if x™ — x™~ x — 1/2 (or x™ = x™Z\ + 1/2 in case the left clock rings) then it remains put. 

2. otherwise it jumps to the right by one and so do all particles x™^ 1 with x™^ 1 = x™ + l/2 for I = 1, 2, ... . 
(in case the left clock rings all particles x™ +l with x™ +l = xT — 1/2 jump to the left for 1 = 1,2,.. .). 
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Figure 1: (a) The initial condition and (b) an example of a point configuration after some time. 

Hence a particle that wants to jump is blocked by particles with lower m-index, but it pushes particles at a 
higher m-index. 

Next we specify the rate of the exponential clocks. For odd m particles jump to the right with rate 
e -1 > and they jump to the left with rate e. For even m the particles jump to the right with rate e and 
to the left with rate e _1 . We are interested in the case where e is small. Hence the particles for odd m 
predominantly try to jump to the right. For even m they mostly jump to the left. However, they are still 
subject to the interlacing condition. Due to this condition, some particles that want to jump to the left (or 
right) are blocked, and even pushed to the right (or left) by particles at a lower level that want to travel 
right (or left). 

If we let time evolve we will see mainly two clouds of particles travelling to the left and to the right 
respectively. See Figure [T] for typical point configuration after some time. For large time we have the 
macroscopic picture as shown in Figure [2] With high probability the particles will be distributed in a 
domain T> contained in the upper half plane. This domain consist of two parts T>\ and T>2- In T>i the 
particles are still densely pact as in the initial configuration (which means the particles did not have the 
chance to jump yet). In T>\ the particles have a density strictly less than one. The boundary of T>\ is a 
smooth curve except for two cusp points. The cusps touch in the limit e \ 0. 

After rescaling the time parameter, the process has a well-defined limit for e | 0. In this limit the particles 
come in pairs. Indeed, by the interlacing condition some particles are blocked and even pushed by a particle 
at one level below that jumps in the reverse direction, thus forming a pair. The pairs are slanted to the right 
in the right half plane and slanted to the left in the left half plane, see also Figure [TJ The process for these 
pairs can be described in the following way: the process decouples in the sense that we have two independent 
processes, one in the upper left quadrant and the other in the upper right quadrant. The process in the right 
quadrant is equivalent to the process where particles can only jump to the right. The process at the left is 
just its reflected version (hence the particle only jump to the left). This process is analyzed in [5] and from 
their results we recover that the limiting domain has a tacnode. 

By standard arguments we can compute the limiting mean density of the particles in all cases. This 
settles the macroscopic behavior for the particles at large time. At the local scale we retrieve the well-known 
universality classes. 

First consider e > 0. If wc zoom in at a point away from the boundary, then we find that the local 
correlations are governed by one of the extensions to the discrete sine process that falls into the class as 
introduced in [3J. If we zoom in at a point at the boundary, but not the cusps points, then we obtain the 
Airy process (see |19j and |10j for a review). The local correlations near the cusps are determined by the 
Pearcey process [21 [8j [9j [TBI [22] ■ Since the proofs of these results follow from standard computations, they 
will be omitted. We do not need these statement for our main results. 

In the case e = we have a decoupled system. In addition to a discrete sine process and the Airy process 
we also obtain the local correlations around the tacnode, which are described by a process that is directly 



related to the GUE minor process (as we will prove). 

The main result of the paper is to derive the process around the tacnode in the transition regime e J. 0. 
The main result is that we obtain a process that has not appeared in the literature (to the best of our 
knowledge). Naturally, it interpolates between the Pearcey process and a process related to the GUE minor 
process. 

In Section 2 we will state our main results and prove them in Section 3. 



2 Statement of results 

We start with some definitions. Let X be a discrete set. A point process on X is a probability measure on 
2 . A point process is completely determined by its correlation functions 

p(X) = Prob{F e 2 X I X C Y}. (2.1) 

A point process is called determinantal, if there exists a kernel K : X x X — > C such that 

p(X) = det[K(x,y)] Xtyex . (2.2) 

For more details on determinantal point processes we refer to [H [TTJ [12l [T4J [T5j [20l [21] . A determinantal 
point process is completely determined by its kernel. 

Now let us return to the evolution on the interlacing particle system as decribed in the introduction. 
By stopping the process at time t, we get a random collection of points on the grid Q. Hence, the Markov 
process at time t defines a point process on Q. In |5j, the authors proved that this is in fact a determinantal 
point process on Q with kernel K given by 

K(x 1 ,m 1 :x 2 ,m 2 ) = ~ Xmi<T " 2 I (1 - £ly )["W2]-[W2] (1 _ £ / w) [( mi +i)/2]-[( m2+ i)/2] w [ Xl ]-[ x ,] ^ 

2?ri J To w 

<W dz — I >— A A-—- - „,.„ -r , (2.3) 



(27ri) 2 / ro / r£ e l e*( Z +^(l -£z)[ m 2/2](! _ e / z )[(m a + l)/2] z [x 2 ] z(u>-z) : 

where [a;] is the largest integer less then x. Here Tq is a contour that encircles the essential singularity but 
not the poles e and e" 1 . The contour T eE -i encircles e and e _1 . Both Tq and r e E -i have anti-clockwise 
orientation and do not intersect each other. Finally, 

_ J 1, if mi < m 2 ,„ .s 

Xmi< " 12 ~ \ 0, otherwise. 1 ' 

To be precise, the variables we use are different from [5 . We choose a symmetric picture since we have 
particles jumping both left and right, whereas the particular model that was analyzed in detail in [S] has 
particles jumping to the right only. Now (|2.3p is obtained by taking the kernel in [51 Cor. 2.26] and substi- 
tuting y = [x] — [(to + l)/2)] and setting the ai's for even values of I to e, and the other ones to e _1 . And 
finally, a conjugation by (— e)[( m2 + 1 )/ 2 ] _ [( mi + 1 )/ 2 )] which does not effect the determinants in the correlation 
functions. 



Our first result is that for large time we obtain the limiting situation as described in the Introduction 
and shown in Figure [2j 

Theorem 2.1. Let H = {z E C | Imi > 0} and F : H -» C given by 

F(z) = r(z + z _1 ) + |log(l + £ 2 -e(z + z- 1 )) - flog*. (2.5) 

Define T>\ by 

£>i = {(£, e » x R+ I 3z £ H F'{z) = 0}. (2.6) 
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(a) x (b) 



Figure 2: The typical shape of the limiting domains T>\ and T>2- The case e > at the left and e = at the 
right. In T>\ the particle are still in the initial configuration. Outside T>\ and T>i there are no particles. 



The boundary dD\ has two cusp points located at 

(0,e + e- 1 ±2). (2.7) 

Set 

t = tL 

m = [fiL] (2.8) 



then the limiting mean density is given by 



1, (£,M)e£> 2 , 



lh^K(x,m,x,rn) = i iargz(£,^i), (£,^)eX>i, (2.9) 
" ' 0, otherwise. 

where z(^,/i) is the unique solution in the upper half plane of the equation F'(z) = 0. 

The proof of this result follows from standard steepest descent analysis on the double integral formula 
for the kernel. Because the proof is standard, we will omit it in this paper. We do not use Theorem 12.11 in 
the sequel. See for a proof of a similar statement in a comparable situation or [16| for an exposition of 
the steepest descent technique on double integral formulas. 

Remark 2.2. The boundary dT>i can be explicitly computed. Indeed, it is clear that 

dT>i = e K x M+ | 3z e R F'(z) = A F"{z) = 0}. (2.10) 

Now for each z G M we have that 

F'{z) = (2 n) 
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Figure 3: A close up picture of the process, for e = 0, around the point where the cusps touch. Note that the 
particle come in pairs that are slanted to the right in the right half and to the left in the left half. Around 
the tacnode, there are forming long vertical strings of particles. 



is system of equations that is linear in £ and [i. So we can easily express (£, /i) as a function of z G R. In 
fact, it is easily checked that for each z G R \ {0, e, £ -1 } the corresponding (£, fi) satisfy £ G R and \x > so 
that (£, /j) is a point on dT>\ . Therefore, the closure of the image of R \ {0, e, e -1 } under this map z i— > (£, /i) 
gives the boundary. The cusp pointing down corresponds to z = — 1 and the cusp pointing up to z — 1. The 
boundary touches the a;-axis at z = e, e^ 1 . 

From Theorem 12.11 it follows that we can achieve the situation of a tacnode in the following way. For 
every fixed e the cusp points on the boundary differ by 4. The location tends to infinity when we take the 
limit e J, 0. By rescaling [i with e the cusp points have a limit as e J. and the gap between the cusp points 
vanishes, resulting in a tacnode. 

From Theorem 12.11 and (|2.3[) we expect to arrive at a process around the tacnode when we scale 



However, it is less clear how to describe the process that arises at the cusp. As shown in Figure [31 there 
will be long vertical strings of particles. In fact, the length of these strings will be of order L and hence the 
density of the particles at this scale will diverge. Therefore we do not obtain a point process in the limit if 
we consider the process on the particles. 

There are several ways of constructing a meaningful process in the limit. Perhaps the most straightforward 
approach is the following: instead of allowing fij to be a free variable, we choose N G N and fix fij G R 
for j = 1,...,N with fij ^ /ife if j ^ k and cut the process at the section rrij = [L 2 (l + /ij/L)]. For 
simplicity, we will also shift the particles on any even m^-section to the left by a half. In this way, we obtain 
a determinantal point process on Z x (1, . . . , TV) for each L with correlation functions 



Prob({particle at (x k , jk) G Z x (1, . . . , N) | k = 1, . . . , n}) = det (K (xk, m jk , x t , m^))^^ (2.13) 




(2.12) 



e = e/L 
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Figure 4: The contours of integration in the kernel /C c 



for all neN. 

Now to obtain the limiting process as L — > oo for the limiting process on Z x (1, . . . , N), it suffices to 
compute the pointwise limit for the kernel K. The limit is given in the next theorem, which is the main 
result of the present paper. 

Theorem 2.3. With t, rrij as in (I2.12[) we have that 

lim K(xi,mi,X2,m 2 ) = K. e (xi, fii,x 2 , H2) ( 2 -14) 

L — >oo 

for (xj , fij) 6ZxR, where 

JC e (x 1 ,Lii,x 2 ,m) = — / e^-^^+^V 1 "* 2 - 1 ^ 

1 r r e ev2(z+±)+4(z+i) 2 w xi dzdw 

+ ttt^t <P <P ; 2 ; 7 7 (2-15) 

{2my J ro 7eus-i e ^(w+^)+^-(w+^) z X2 z \ w ~ z ) 

The contours of integration and their orientation are as indicated in Figure^ More precisely, Tq is a contour 
encircling the origin with counter clockwise orientation. The contour £ is a contour connecting — ioo to ioo 
that does not intersect To and stays in the right-half plane. 

To the best of our knowledge, the kernel K e has not appeared in the existing literature yet. 

It is not difficult to show that after inserting the new parameters given in (|2.12p in the integrands in 
(|2.3[) and taking the pointwise limit as L — > oo, one obtains the integrands as given in the right-hand side of 
(|2.15p . However, there is an important technical issue that needs to be taken care of. Note that the integrand 
in the double integral contains poles at s and e -1 , and also an essential singularity at 0. By taking the limit, 
the pole approaches the essential singularity at the origin, which complicates the contour deformation in the 
analysis. 

A different way of creating a point process is the following. Instead of considering the location of the 
particles, one could consider the statistics of the upper endpoints of the vertical strings of particles. We will 
restrict our process only to the odd m-sections. Then (x, m) is defined to be a upper endpoint if there is 
a particle at (x,m) but there is no particle at (x,m + 2). The upper endpoints form a point process with 
correlation functions 

Pn{{x\, mi), . . . , (xn, mjv)) = Prob (particle at (Xj, rrij) and no particle at (xj, rrij + 2) | j = 1, . . . , N) 

(2-16) 

for (Xj,mj) € Z X (2N - 1) and NeN. 



G 



Theorem 2.4. With t = eL and rrij the closest odd integer to L 2 (l + fij/L), we have that 

lira L N p N {(xi,m{), . . . , (x N ,m N )) = e N det (IC e (xi - l,fii,Xj,Hj) + K e (xi + 1, fa, Xj, /Zj-)) . . (2.17) 
for (x j ,» j )GZxMforj=l,...,N. 

Remark 2.5. Another argument for the fact that the lengths of the vertical strings of particles must be of 
order L, is that one can prove that the process restricted to two horizontal sections that are close, is just 
constructed of two copies of the process restricted to one of these horizontal sections. To be precise, let 
/ii £ R and take mi = [L 2 (l + pi/L)}. For each L let m 2 be such that m\ — m 2 — o(L) as L — > oo. For 
simplicity, we assume that m 2 > m,\. Then it is not difficult to prove that 

lim K(xi, mi, x 2l mi) — IC £ (xi, fii, x 2 , pi), (2.18) 

L — >oo 

lim K(xi,m 2 ,x 2 ,m 2 ) — JC e (xi, pi,x 2 , pi), (2.19) 

L — >OG 

lim K(xi,m 2 ,x 2 ,mi) = IC e (xi, pi,x 2 , pi), (2.20) 

L — >OG 

lim K(x 1 ,mi,x 2 ,m 2 ) = -6 Xl X2 + K e (x%, fix,X2, fa), (2-21) 

L — >oo 

for all x\,X2 E Z. Hence, for a^i , . . . , x\ , t/i, . . . , we have 

lim Prob(particles at (a^mi), . . . , (x/,mi), (j/i,ma), ■ ■ ■ , {yk->rn 2 )) 

L— >oo 

= det( (^(^'Mi^Mi))^! (^ e (^,Mi,%,Mi) -^.x 3 )Jf lj=1 \ ^ 22 ^ 
\j£%i,Mi,^pMi))ti,j=i (^ e (j/i»Mx,%-)Mi))* 3 - = i • / 

This implies that (in the limit L — » oo) the process on the line m2 is just a copy of the process on the line 
mi. 

We will now derive some properties of the kernel /C e . The following proposition shows the symmetry in 
the kernel. 

Proposition 2.6. We have that 

1. £ 6 (— sgi,/Zi, -x 2 ,p 2 ) = K, e (xi - l,fa,x 2 - l,/i 2 )- 

& (-l)* 1- * 3 *^!, -/ii, X 2 , -/ia) = fyxi^iMza./w) ~ ^(zi, Mi, ^2, M2) 
/or all (xj,/J,j) 6ZxM. 

The first property shows that our point process is invariant with respect to the transform m-l-i. 

To interpret the second symmetry property in this proposition, we note that if V is a determinantal 
point process on a discrete set X with kernel K, then we have that 1 — K is the kernel of the determinantal 
point process V' defined by V'(X) — V(X \ X), for X C X. This process is sometimes referred to as the 
complementary process. It is obtained by replacing particles with holes and vice versa. For more details on 
this particle hole involution we refer to the appendix of [BJ. 

In the final results of this paper we investigate the limiting behavior of the kernel as e J. and e — > 00. 
We start with the first case. 

Theorem 2.7. For x\ : x 2 < we have that 



lime* 1 X2 JC e (x 1 ,fii,x 2 ,fi 2 ) = -x^k^Xx!^^- IJ-iY 2 Xl 

ej.0 



1 f f e ^ z+ 2 z w Xl dzdw 
+ (2^ t 7 E e ^-+^ 2 w-z' ( 3) 
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Figure 5: From the GUE minor process to the point process with kernel given in (|2.24[) . The open circles 
represent the point {x, y} x ) from I = 0, . . . , x and x = 0, 1, . . .. In the left picture we draw the vertical lines. 
The dotted lines are only auxiliary. In the right picture, we draw the lines associated to the choice of the \ij . 
The solid circles are the intersection point of the dashed horizontal and solid vertical line. The solid circles 
describe the process with kernel (|2.24l) . 



and for X\,X2 > 0. 



lime* 2 Xl tC e (x 1 ,ni,x 2 ,p.2) = -X l n< t i 2 Xx 2 <x 1 (^2 - Mi)* 2 

<40 



1 f f C ^ Z+ 2 Z z X2 dzdw 

(2.24) 



(27ri) 2 J Tq J s e viw+h w2 w Xl+1 w - z 
Moreover, the limit of the kernel vanishes if neither of the two conditions on x± and Xi are satisfied. 

The one dimensional version \i\ = \ii of this kernel has appeared before in the literature, see [TJ. By 
expanding the term (w— z)^ 1 we can express the double integral as a sum of products of Hermite polynomials. 

By inserting (|2.24|) (or (|2.23|) 1 in (|2.17[) we find the limiting process for the endpoints in the right (or left) 
half plane. It turns out that in each of the half planes one of the kernels in the determinant tends to zero 
as e I 0. In fact, the limiting kernel is the kernel corresponding to the GUE minor process, see for example 
[T51 117j . Indeed, in the case x\,x 2 > we rewrite (12.171) as 

lim L N pM^xx^mx), (x N ,m N )) = dot ( e *i+i-*< fc e ( Xl - 1, fx^Xj, Hj) + e 2 e l>i ~' l ~ xl K e {x i + 1, /x i( x jt ^j)) 

(2.25) 

and then by (|2~24]) we find 

lim lim L N p N ((xi, mi), . . . , (x N , m N )) = det (K GVE (xi, pi, x 3 ■, Hj))f 4 , (2.26) 

elO L—*oc ' J 



where 



3 ^ z +i z Z X2 dzdw 



The case X\,x 2 < is similar. By comparing with [T31 Def. 1.2] we see that the kernel Kque describes the 
GUE minor process. 
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Figure 6: The contours of integration for the Pearcey kernel. 



Remark 2.8. Based upon the relation between the construction of the two processes just above and below 
Theorem l2.3[ we see how the process defined by the kernel at the right-hand side of (|2.24j) can be constructed 
out of the GUE minor process explicitly: pick a point configuration random from the GUE minor process 
and denote the points by (x,y l x ) £Nx R. At each vertical x section we draw x vertical line segments. Each 
segment has y l x as a lower endpoint and y l ~\ as the upper endpoint. Here we define y® = oo. See also Figure 
We next fix N different real numbers fij € R. Then we define a point process in N x {1, . . . , N} by 

Probjparticle at (xi,ji), i = 1, . . . k)} — Prob{Each (xj, is on a line segment} (2.28) 

The conclusion is that the new process is in fact a determinantal point process with kernel as in (|2.24|) . 



The second situation that can be obtained is that for e — » oo, in which the cusps should separate. Hence 
we expect to obtain the Pearcey process when we take the simultaneous limit fij — > — oo (or fij —* +oo). 
The following Theorem states that we indeed find the Pearcey process at the lowest of the two cusps. 

Theorem 2.9. Set 

e = M 

(ij = -M(l - Vj/2M) (2.29) 
Then 

J im mi/2 fC e {xx,^,x 2 ,^) = / dw 

M^oo M l l Z 27T1 ./_:_, 



1 r io ° f e i zi +^ z2 -^ z dwdz 



(2.30) 



(27ri) 2 J_ ioo J c e 2 w4 +^i w2 -ii w w - z 
The contour C conists of four rays, from ±e m / 4 oo to and from to ±e 3,n / 4 oo (see also Figure^. 

The kernel given (|2.30p is known in the literature as the extended Pearcey kernel that describes the 
Pearcey process, see [2 HI [9J [HI [22] for more details. 

To conclude this section, note that if we combine Theorem 12.91 with the second symmetry property in 
Proposition 12.61 then we see that the complementary process near the top cusp locally converges to the 
Pearcey process, as expected. 



3 Proofs 

In this section we prove our results. 
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3.1 Proof of Theorem EO 

Let us first introduce some notation. Write the kernel in (12.31) as 



TSl \ Xrn 1 <m 2 I Gt miXl [Z,)Az ^ I I Gf ,tni ,xi \W) 1 , , /„ 

K(xi,m 1 ,x 2 ,m 2 = — — f — — —r^— + 77r^f f n hw 7 dzdw, 3.1 

27T1 Jr G t,m 2 ,x 2 (z) 2 (27T1J J JtqJt^! G t,m 2 ,x 2 {z) Z(W - Z) 

where 

Gt lWIj x(ifi) = e t{w+1/w) (l - ew)^' 2 Hl - e/wf mi+1) /^w^ . (3.2) 

Let us for a moment ignore the contours of integration and take the limit L — > 00 for the integrand. It is 
not difficult to see that we have the following pointwise limit 

Q. ( w ) e e t i 2 ( Z +l/z)+e 2 (z+l/z) 2 [xi] 

L^o G t>m2:X2 ( Z ) ~ e ^2(w+l/w)+eHw+l/ w r «[x a ] ■ ^- 6 > 

for z, w S C \ {0}. Hence the difficulty in the proof lies in the fact that we have to control the contours of 
integration (which clearly depend on e and hence L) while taking the limit. 
As a first step we prove the following result. 

Lemma 3.1. With G as in (13.21) and K as in (12.31) we have that 



K {x 1 ,m 1 ,x 2 ,m 2 = t-^— * dz — — + —-7^ * f — r (3.4) 

Here T £ and T s -i are two contours encircling the poles e and e^ 1 respectively, but no other poles. The 
contour To, £ is a contour encircling the origin and the contour T £ but not T £ -i. The contours are taken so 
that they do not intersect and have counter clockwise orientation. See also Figure [?| 

Proof. First split the contour T £ e -i into two small contours T £ around e and T E -i around e~ 1 . Then 



Xm 1 <m 2 f Gt miXl {z)6.Z Iff Gt m i,xi(w) 1 , , 

K(xi,m 1 ,x 2 ,m 2 ) = — — * — - — '— r -— + jtt-^ f f tt — —TT~( \ dzdw 

2tti Jr Gt, m2 .x 2 \z) z [^ly J To J Ts & t , m2)X2 (z) z(w - z) 

1 II G t ,mx,xi{w) 1 



I G t , muXl {w) 1 
f 7; ■ — TT~( \ dzdw ( 3 - 5 ) 

Jr Gt.rn 2 .x 2 (Z) Z(W - Z) 



(2m) 2 J ro J Ve l G t . m2 ,x 2 {z) z(w - z) 

Now in the first double integral in the right-hand side we deform To so that it encircles T e . This deformation 
will be denoted by Tq^ £ . Note that we now pick up a residue at w = z. Hence by deforming we obtain an 
extra single integral over the contour T e 

tsi \ \m\<im 2 f Gf ffi l X ,(zj dz 1 f Gt ?m ,x\ (z) dz 

K(xi,mi,x 2 , m 2 ) = — - 



2-ki J ro G t ,m 2 ,x 2 (z) z 2ni J Fs Gt, m2 ,x 2 (z) z 
Gt — H 1 dzd w + T ^{ I G ; r ^] 1 dzdw (3.6) 



(27ri) 2 J To £ J Te G t ,m 2 ,x 2 (z) z(w - z) ' (27ri) 2 J To J v ^ G t , m2 ,x 2 (z) z(w - z) 

Now the extra single integral has the same integrand as the first single integral. The integration is now over 
a contour encircling the pole at e (and no other pole). However, this pole is not present in the case mi < m 2 
and then the integral vanishes. Therefore we can glue the integrals over To and r e together and obtain 

Ti r, \ Xmi<m 2 f G t . mi .xi ( z ) dz 1 / / Gt mixi (w) 1 

K(xi,m 1 ,x 2 ,m 2 ) = — — * — --— + f f — --— -dzdw 

2tti Jt , c Gt,m 2 ,x 2 (z) z (2tt\Y J T ^ £ J Te G t ,m 2 ,x 2 {z) z(w - z) 

1 IS Gt,mi,xi(w) 



f n ( \ ~< T dzdw ( 3J ) 

Jr _i Gt,m 2 ,x 2 {z) z(w - z) 



(2?ri) 2 J To J Ve l G t , m2 ,x 2 {z) z(w - z) 

Finally, note that the integrand has no pole at w — s so in the second double integral we can safely 
deform Tq to be ro, £ . This proves the claim. □ 
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Figure 7: Deforming the contours as in Lemma 13. II 

Now the proof of Theorem 12.31 follows by simply taking the limit in the integrands and correctly choosing 
the contours T e and r e -i. 

Proof of Theorem \2.3[ Since e will be small, we take the contour Tq i£ to be some fixed contour encircling 
the origin independent of e, say the unit circle for simplicity. Due to the behavior of Gt t m 2 ,x 2 ( z ) f° r z ~* 00 
we can deform T £ -i to any simple contour £ that connects 00 to 00 crossing the real axis once at a point 
between 1 and e _1 and contained in the sector 

{zeC \ -it/2 + 5 < argz < tt/2 - 8}. (3.8) 

for some fixed 8 > 0. In fact we choose S to be a contour so that it eventually falls inside the sector 

{zeC \ -tt/2 + <5 < argz < -tt/4 - 8, vr/4 + 8 < argz < tt/2 - 8}, (3.9) 

for some 8 > 0. 

The behavior of Gt l1Tl2>X2 (z) for z — > is to a large extent similar to the behavior near 00, which can be 
seen by performing the transform z 1— > 1/z. We can (and do) deform the contour T e to be E _1 . Note that in 
this way, we have deformed the contours, so that they go through the essential singularities of the integrand. 

Now that we have chosen the contours, that are clearly independent of e, we can compute the limits of 
the integrand, which is given in (|3.3[) . Due to the choice of the sector in which £ eventually falls inside, the 
convergence is uniform in z e E U E_i and w G IV This proves the statement. □ 



3.2 Proof of Theorem Q 

Proof. We start by noting that by expanding a factor (1 — e/w)(l + ew) in the integrals for the kernel K in 
(|2.3|) we get for odd n, m 

K(x, n + 2 7 y, m) = (1 + e 2 )K(x, n, y, m) — e(K{x + 1, n, y, m) + K{x — 1, n, y, m)) 

+ %^l(l-Eio)(l-£/!i;K- ! '- 1 d(D. (3.10) 

27T J 

There is an additional single integral because of the fact that £ n +2<m = Cn<m — <5n,m-i — 8 n , m -2- Now since 
n, m are both odd 8 n ,m-\ vanishes trivially and only 8 n , m -2 remains. 
Similarly, 

K(x, n, y, m — 2) = (1 + e 2 )K(x, n, y, m) — e(K(x, n, y + 1, m) + K(x, n, y — 1, to)) 

+ On^zl Li _ ew )(i _ e/w^-v-ifao. (3.11) 
27r 7 
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Therefore 

— K(xi, riii + 2, Xj, rrij) + K(xi, m t , Xj, rrij) — e\K{xi + 1, rrii, Xj, rrij) + K(xi — 1, mi, Xj,rrij) 

) + 0(e 2 ) (3.12) 

Now note that if we have rrii ^ rrij then m; — rrij is of order L, so that we can ignore the term £ mi , m --2i so 
that 

— K(xi,rrii + 2,Xj,rrtj) + K(xi,mi,Xj,rrij) = e( y K(xi + l,mi,Xj,mj) + K(xi — 1, rrii, Xj, rrij) j (3.13) 
and 

— K(xi,mi + 2 1 Xj, rrij + 2) + K (xi, rrii, Xj,nrij + 2) = e(^K (xi + l 1 mi,Xj, rrij + 2) + K{xi — 1, m», xj , m 3 - + 2) 

) + C?(e 2 ) (3.14) 

By subtracting (|3.14[) from (|3 . 13[) and inserting (|3.11[) we obtain 

—K(xi, rrii+2, Xj,mj)+K(xi,m l , Xj,rrij)+K(xi, rrii+2, Xj,rrij+2)-K(xi, mi,Xj,mj+2) = 5 mi ,m j 5 XuXj +0{e ) 

(3.15) 

By the complementation principle (see [6]) 

MO^mx), . . . , (s^m*)) = dot (^ (xi)m . ja; ./ TO . + 2) J - Kix^m, + 2,x J ,m- + 2) ) ^ 

By adding the first column to the second and afterward subtracting the second row from the first we obtain 
that the determinant equals 



det 



( K(xi, rrii, Xj, rrij) — K(xi, rrii, Xj,rrij + 2) —K(xi, rrii + 2, Xj, rrij) + K(Xi, rrii, Xj,mj) \ 

-5( Xi , mi ),( Xj , mj ) + K(x l ,m l + 2,Xj,rrij + 2) - K{x il m il Xj,mj + 2) 



\ K(Xi,m,i,Xj,m,j +2) ^(a i ,m i ),(x ;j ,m :f ) - K(xi, m l + 2, xj, rrij + 2) + K(xi, rrii, xj, rrij + 2) J 

(3.17) 

Now using, (|3.1ip in the upper left block, (|3.15|) in the upper right block, and (|3.14[) and (|3 . 1 1 [) in the lower 
right block, this determinant has the form 

p N ((x 1 ,m 1 ),...,(x N ,m N ))=det^ m £{K{Xi + ljmiiX3jmj) + K{x ^ ljmuXjjmj)) ), (3-18) 
from which the proposition easily follows. □ 



3.3 Proof of Proposition 12.61 

Proof. 1. The statement easily follows by the transform z i— > 1/z and w i— > 1/w in the integral representation 

Of K e (—Xi,fXl,-X2,H2). 

2. The second property follows by the transformation z i— » — z and w >—>■ —w in the integrals and a 
deformation of S. In the definition of the kernel if e we take S on the right of r . However, we can also take 
S to be on the left at the cost of an extra integral as shown in Figure [5] The second double integral at the 
right is easy to compute since the integral over z encircles the pole z — w only and hence can be computed 
by the Residue Theorem. The result of this is that we can rewrite the kernel in the following way 

K e (xx,fXi,x 2 ,fJ,2) = S^ ulluX3ilia ) +£m2</*i / ~f / (3-19) 

Jr Jr J-(sus- 1 ) 

Now by inserting fj,j = —fij and the transformation w i— > —w and z i— > — z we arrive at the statement. □ 
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Figure 8: In the definition of JC e we can switch the contour £ to — £ to obtain (|3 . 19[> . 



3.4 Proof of Theorem [27T1 

Proof. By the first symmetry property in Proposition 1 2 . 61 it suffices to consider the case x-i < only. Using 
the transform z i— > z/e and i» h ui/e we obtain 



1 /" /" e^+^+^+^'ur 11 dzdw 



r 



(27ri) 2 / ro y suS -i e Mi(t"+^r)+5(' i, +i7r) 2 z 2;2 - z) 

Setting e = at the right-hand side gives 



(3.20) 



\ime Xl ' X2 ^(x 1 , f i 1 ,x 2 ^ 2 ) = -x Ml<M2 / e^-^^w^-^-Mti; 



1 /" /" e M2z+2^ 2 u , a: i dzdw 



r 

+ (27ri) 2 % a 7 sus -i e w w +>V 2 z(w - z) ' ^ 3 ' 21 ^ 

The single integral can be easily computed. As for the double integral, we note that since x 2 < we have 
that the integral over E" 1 vanishes and hence 

i r r c^z+^z^x, dzdw 

(27ri) 2 J To 7 S e W™+§«>V 2 z(w - z)' 

The single integral can easily be computed. This proves (|2.23p . 

Finally, note that if x\ > then both integrals vanish since there is no pole for w = 0. □ 

3.5 Proof of Theorem EH 

Proof. The double integral in (|2.15|) in the new parameters given by 12 . 29[) reads 

e ^(z-l) i /z 2 +M0 2 (z+l/z)-M 1 ' 2 Z 1 lna 1 

dw I dz — = ; (3.23) 



(27Ti) 2 Jr Q JsuS- 1 e ^-(w-l)*/z 2 +M& 2 (w+l/w)-M 1 / 2 & In a z(w - z) 
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Figure 9: The line {w | Im((w — l) 4 /w 2 ) — 0} which are the contours of steepest descent /ascent leaving 
from w = 1. 



For large M the main contribution comes from the terms 

^ ^- and i 5^-. (3.24) 

z w 

This term has a critical point at z = 1 (resp. tu = 1) of order three. Therefore, there are four paths of 
steepest descent leaving from z = 1 and four paths of steepest ascend, as shown in Figure [9] In fact, the 
paths of steepest ascent leaving from w = 1 are the unit circle and the real line. We now deform r to be 
the unit circle. The contour £ is deformed so that it passes through z = 1 and follows the path of steep 
descent outside the unit circle. Then locally around z = 1 and w = 1 the contours Tq, S and E _1 can be 
deformed to the contours for the Pearcey kernel as shown in Figure [S] 

By standard steepest descent arguments one can now show that the dominant contribution comes from 
a neighborhood around z = 1 and w = 1. More precisely, after introducing the local variables 



w = 1 + w/M 1 / 2 
z = l + z/M 1 ' 2 



(3.25) 



it is not difficult to prove that 



If f e M?-(z-l) 4 /z 2 + Mis 2 {z + l/ Z )-M 1/2 t2ln Z 

dw / dz 



(27Tl) 2 Jr Q JSUS- 1 e ^-(w-l) i /z 2 + ML> 1 (w + l/w)-M 1 / 2 ti 1 \nz z(w - z) 

dz / dw 



(27ri) 2 j-ioo Jc e 2* 4 +^n5 2 -?i«' w 



(^(A/- 1 / 2 )) , (3.26) 



as M — » oo. 

Now consider the single integral in (|2.15[) . which in the new parameters (|2.29p reads 

— I dw e^-^+V^w™ 1 ' 2 ^-^- 1 (3.27) 
2tti J To 

The dominant term in this integral is 

e M(v 2 -vi){w+l/w) _ (3.28) 
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A simple computation shows that w + 1/w has two saddle points w = ±1, both of order one. Since V2 > V\ 
(otherwise the single integral is not present) we have that w = +1 is the dominant saddle point. This means 
that, the main contribution comes from the part of the contour close to w = 1. Hence, if we introduce the 
new local variable 

w = l + w/M 1/2 , (3.29) 
then by standard arguments one can prove that 

— I dw eMMK"^)/''^)-! 

ul/2 2M(i/!-ki) rioo 

= — / dw e te-»i)*'-fr-b)* ( 1 + OiM- 1 ' 2 )) (3.30) 

2nl J-ioo V ' 

By inserting Q3.30P and (|3.26|) in ()2 . 1 5[) and taking the limit M — ► oo we obtain (|2.30p . This proves Theorem 

HI □ 
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